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SUMMARY 


A new trigonometric approach to the finite difference calculus was applied 
to the problem of beam buckling as represented by virtual work and equilibrium 
equations. The trigonometric functions were varied by adjusting a wavelength 
parameter in the approximating Fourier series. Values of the critical force 
obtained from the modified approach for beams with a variety of boundary condi- 
tions were compared to results using the conventional finite difference method. 
The trigonometric approach produced significantly more accurate approximations 
for the critical force than the conventional approach for a relatively wide 
range in values of the wavelength parameter; and the optimizing value of the 
wavelength parameter corresponded to the half-wavelength of the buckled mode 
shape. Thus, selection of the wavelength parameter is a simple process if the 
half-wavelength is known. Methods for selecting this parameter in the more 
general case are also presented. It was found from a modal analysis that the 
most accurate solutions are obtained when the approximating function closely 
represents the actual displacement function and matches the actual boundary 
conditions. It is more difficult to select a satisfactory value of the wave- 
length parameter for the equilibrium equation which makes the virtual work 
equation more attractive for practical applications. A comparison of the vir- 
tual work and the Galerkln approaches identified marked similarities between 
the two methods. 


INTRODUCTION 


Numerical analysis has played an important part in furthering the under- 
standing of structures over the past decades. Finite element and difference 
are just two techniques which may be considered part of the overall numerical 
approach o The similarities between the two have come to light recently as 
they are applied to energy equations, and therefore further research into 
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finite difference methods, with respect to structures, has been pursued. Stein 
and Housner (ref. 1) initiated work into a trigonometric approach to finite 
difference as applied to plate instability which appeared very appealing be- 
cause of its convergence characteristics (ref. 2). Thus, the authors investi- 
gated this relatively new trigonometric approach and applied it to problems of 
beam buckling. Incorporating both the virtual work and equilibrium equations by 
extending the method's fundamental mathematical concepts. In order to accom- 
plish this, a trigonometric function was varied by adjusting a wavelength par- 
ameter in an approximating Fourier series. Values of the critical force ob- 
tained from the modified approach for beams with a variety of boundary condi- 
tions were, for the first time, at least to the authors' knowledge, compared to 
results using the conventional finite difference method. The trigonometric 
approach produced significantly more accurate approximations for the critical 
force than the conventional approach for a relatively wide range in values of 
the wavelength parameter; and the optimizing value of the wavelength parameter 
corresponded to the half -wavelength of the buckled mode shape. It was found 
from a modal analysis that the most accurate solutions are obtained when the 
approximating function closely represents the actual displacement function. It 
is more difficult to select a satisfactory value of the wavelength parameter 
for the equilibrium equations which makes the virtual work equation more 
attractive for practical applications. The buckled mode shape (or eigenfunc- 
tion) is predicted with high accuracy regardless of the value of the wavelength 
parameter. A comparison of the virtual work and the Galerkin approaches iden- 
tified marked similarities between the two methods. 


MATHEMATICAL EQUATIONS AND NUMERICAL APPROXIMATIONS 


Virtual Work Equation . The virtual work principle can be stated mathe- 
matically as 

6W = 6U (1) 

e 

(We = external work; U = strain energy) 


For a one-dimensional beam with an axial force P, this relationship reduces to 
the following equation (ref. 3) ; 
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where v represents the vertical displacements during buckling. The derivatives 
in equation (2) can be replaced by trigonometric finite difference approxima- 
tions : 
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where h = mesh spacing and X = buckle wavelength parameter. (Refer to the 
Appendix for complete derivation of equations (3) and (4)). Note that the vir- 
tual work equation has been evaluated using half-station approximations for the 
first derivative and full-station approximations for the second derivative. 

The trigonometric finite difference expressions represented by equations (3) 
and (4) reduce to the conventional polynomial expressions as X approaches in- 
finity since 


lim ■rr _ ^ (5) 
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Integration in equation (2) is performed using the trapezoid rule. The virtual 
displacements in the resulting equation can be ordered to produce 

N 

I f.(v,P)5v. = 0 (6) 

i=l ^ 


For this equality to hold, the coefficients of the individual 6v . must be equal 
to zero which leads to an eigenvalue problem of the form 

N 

I g..(P)v =0 j = 1,2,---N (7) 

i=l 1 

from which the critical force can be calculated. 


Equil i brium Differential Equation. The elastic form of the differential 
equation describing the slightly bent equilibrium configuration of an initially 
flat beam can be stated as (see ref. 4) 


d^v P d^v _ 

dx^ El dx^ 


( 8 ) 


The trigonometric finite difference approximation for the fourth derivative was 
derived incorporating a five term Fourier series with the following result: 


where 
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= -2 cos 6+2 
= 2 cos(29) - 2 
A^ = -4 cos(20) + 4 cos 0 


( 12 ) 

(13) 

(14) 


A^ = 16 sin^0 - 48 cos^0 + 32 cos^0 + 48 cos^0 - 32 cos^(0) (15) 

Equation (8) can b.e evaluated at each of the N grid points placed along the 
length of the beam, resulting in a set of equations with the following form: 


— ;r(v. „ - 4v. , + 6v. - 4v + v. ,„) 
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(16) 


where A represents the combined coefficient resulting from the evaluation of 
equation (9) . The displacements in this set of equations can be ordered to 
produce an eigenvalue problem in the form of equation (7) . 


Virtual Work and Galerkin Approaches . To complete the discussion of fi- 
nite differences, it is instructive to compare the virtual work approach, as 
applied in this paper, to the Galerkin approach. Both methods are based on the 
principle of minimum potential energy. The Galerkin approach approximates the 
displacement function by the following series expansion: 

M 

v(x) = y a.(f). (x) (17) 

.1 11 
1=1 

where a^ are undetermined coefficients and (|) . represent continuous functions. 

By substituting this approximating function in the potential energy expression 
and performing a sequence of variational operations, the following system of 
equations is obtained: 

M 

y f (P)a =0 i = 1,2 M (18) 

1=1 3 

where f^j denotes a functional relationship in terms of the external force. 

The concept underlying the Galerkin approach is based on the fact that the 
error in the approximation of equation (17) is minimized for any value of M if 
the a^^ are chosen such that equations (18) are simultaneously satisfied (ref. 

5) . The virtual work approach incorporating finite difference approximations 
can be described in a similar manner. The Fourier coefficients, T^, are im- 
plicitly selected (through the computational properties of the virtual work 
algorithm) such that the error due to the approximation of v(x) is minimized. 
This hypothesis was substantiated numerically by an analysis of the Fourier 
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series approximating function. 


CONCLUSIONS AND RESULTS 


This paper compares the results of using the trigonometric and convention- 
al approaches to the finite difference calculus in order to solve the equili- 
brium and virtual work equations. A wide range of boundary conditions were 
Investigated. Various values of A were used in the trigonometric approach to 
determine the optimum value as well as to determine the range over which the 
trigonometric approach gives more accurate approximations than the conventional 
approach. In addition, an in-depth search was conducted to provide plausible 
explanations for the superiority of one method over the other and one value of 
X over other values. Finally, the effect of decreasing the number of grid 
points and the use of full- station approximations in the virtual work equation 
were investigated. 

The virtual work method was found to be an efficient and simple approach 
and provided excellent results for both the trigonometric and conventional 
techniques with as few as five grid points. As predicted from theory, compu- 
tational data revealed that the magnitude of error in computing varies 
directly with the square of the grid size. The variable input parameter, X, 
has the effect of adjusting the wavelength of the Fourier series approximating 
function, and the optimum value of A corresponds to the half-wavelength of the 
buckled mode shape for each boundary condition. There is a range in the val- 
ues of A for which the trigonometric approach is superior. This range extends 
from approximately 25% below the optimum value to infinity. Thus, a large 
value of A is guaranteed to provide more accurate results than the conventional 
approach. Of course, if A is chosen to be too large, the error from the con- 
ventional and trigonometric techniques approaches the same value; and the benefit 
of using the trigonometric technique is lost. It can be shown that A=1,5L 
produces satisfactory results for all boundary conditions. An illustration of 
this can be observed in figures 1 and 2 which depict the error in the calcu- 
lation of P^j- for plnned-pinned and free-guided beams. A potential explanation 
for the superiority of one method over the other was found from an investiga- 
tion of the Fourier and Taylor series approximating functions. The A value 
which yields the closest series approximation to the theoretical displacement 
function corresponds to the optimum value of A. That is, the key to calculat- 
ing an accurate estimate of the critical force is to supply an approximating 
function which very closely reproduces the buckled mode shape as well as to 
satisfy the geometric and force boundary condition. 

The similarities between the virtual work technique as employed in this 
paper and the Galerkin approach have been explored. In both cases, equilibrium 
expressions are used to derive potential energy relationships; and the dis- 
placement functions are approximated by series expansions. In addition, there 
is a strong relationship between the resulting sets of equations developed by 
the two methods. Both methods attempt to minimize the error in approximating 
the displacement function. When the approximating function is altered such 
that this minimized error is larger, the error in the computed critical force 
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will increase proportionally. This concept was demonstrated by the use of the 
full-station finite difference approximation for the first derivative. The 
decreased accuracy in the approximation of v(x) caused a significant increase 
in the error of the calculated value of the critical force. 

The equilibrium approach was also found to be efficient, and excellent re- 
sults were obtained using the trigonometric technique. Figures 3 and 4 depict 
the results of calculating the critical force for pinned-pinned and free- 
guided beams. However, it is more difficult to select an effective value of X 
using this approach than was found to be true for the virtual work approach. 
There are two optimum values of X due to the fact that five terms were used in 
the Fourier series for the derivation of the fourth derivative. The first 
value corresponds to the half-wavelength of the buckled mode shape. A precise 
estimate of the buckled wavelength is required in this case, however, since 
there is little margin for error. The range around this optimum value for 
which the trigonometric approach is superior to the conventional approach is 
very small, and the error builds rapidly as estimates of the optimum value 
worsen. There is a large range around the second optimum value for which the 
trigonometric approach is superior. This range extends from approximately 27% 
below the optimum value to infinity. This provides a comfortable margin of 
error for selecting X. The problem is that there is no known physical parame- 
ter from which this second optimal value can be estimated other than the rea- 
lization that each term of the series is superimposed upon each other, yielding 
a wavelength equal to the theoretical value at this particular X. It appears 
from the available data that a value which is 2.75 times the half- wavelength 
provides a reasonably close estimate in most cases, but specific boundary con- 
ditions vary considerably from this figure. Despite the uncertainty, it is 
much safer to attempt an estimate of the second optimal value of X due to the 
larger error margin. An attempt to use the first optimal value is probably 
unwise unless the buckled mode shape is known a priori with reasonable accuracy. 
It can be shown that X=3.75L provides more accurate results than the conven- 
tional approach for all boundary conditions. 

In comparing the results of the virtual work and equilibrium approaches, 
many similarities were noticed despite the major conceptual differences in the 
derivation of these methods. The interpretation of the wavelength parameter, 

X, is the same in both cases as already discussed. In addition, the virtual 
work and equilibrium methods give the same value for Pj,j- when conventional 
finite difference expressions are used. The two methods do not give the same 
result when trigonometric expressions are used due to the presence of the two 
additional Fourier series terms in the equilibrium equation. Several major 
differences were also noted in the two methods. For example, it is more diffi- 
cult to predict the optimum value of X for the equilibrium approach. Addition- 
ally, it was found that an error in the estimate of X produces a larger error 
in the computed value of P^j. for the equilibrium method than for the virtual 
work method. For these reasons, the virtual work method is recommended for 
general use over the equilibrium method. The trigonometric approach to the 
finite difference calculus is recommended over the conventional approach, par- 
ticularly in those cases when the shape of the displacement function is known 
within rather broad tolerances. 
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APPENDIX 


Trigonometric finite difference approximations are derived in a manner 
similar to the conventional expressions with the exception that the following 
form of the Fourier series is used rather than the Taylor series: 


v(x) 


= T, 


+ T2sin 


Tr(x - Xq) 
"X 


+ T^cos 



(19) 


The derivative of equation (19) evaluated at the reference point, x^, is given 
by 


and 


= ^2 X 

^2 = v'(x„) A 


Evaluation of equation (19) at + h/2 and x^ - h/2 results in 
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( 20 ) 

( 21 ) 

( 22 ) 

(23) 


Subtract equation (23) from equation (22) to obtain 

. TTh 

''-H, - ''JJ - 2J- 


(24) 
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If equation (21) is substituted in equation (24) and the terms rearranged, the 
following expression is obtained 


where 


= T 
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(V. 


4^2 

^ 2Asin(||) 

7T 



(25) 

(26) 


The trigonometric finite difference approximation for the second deriva- 
tive can be obtained in a similar manner. Equation (19) can be evaluated at 
X = Xg + h and x = x^ - h to provide 

irh Trh , ^ , 

v^j^ = + T2sin — + T^cos -y (27) 
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By adding equations (27) and (28) and subtracting two times equation (19) 

obt 

irh 


evaluated at x = x^, the following expression is obtained: 


v+i - 2v^ + v_j = 2T3 (cos ^ - 1) 


(29) 


The second derivative of equation (19) with respect to x evaluated at x = x 


IS 


Solving for yields 


V = -T — cos 

° 3 .2 X 


"3 = - O \ 


(30) 


(31) 


If equation (31) is substituted in equation (29) and the resulting expression 
solved for v^, the following expression is obtained: 


V = 


° 4x2sin2(^) 

2X 
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or 
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Note that the trigonometric and conventional finite difference expressions 
are similar with the mesh spacing, h, simply replaced by h. 
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FIG. 


FIG, 2 
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ACCURACY COMPARISONS FOR USING THE VIRTUAL WORK EQUATION. 
PINNED-PINNED BOUNDARIES. 



ACCURACY COMPARISONS FOR P^p USING THE VIRTUAL WORK EQUATION. 
FREE-GUIDED BOUNDARIES. 
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FIG, 3 ACCURACY COMPARISONS FOR USING THE EQUILIBRIUM EQUATION. 
PINNED-PINNED BOUNDARIES. 



FIG,- 4 ACCURACY COMPARISONS FOR P^r USING THE EQUILIBRIUM EQUATION. 
FREE-GUIDED BOUNDARIES. 
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